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The stability of the brick wall model is analyzed in a rotating background. 



£> ' It is shown that in the Kerr background without horizon but with an inner 



boundary a scalar field has complex-frequency modes and that, however, the 



(^ , imaginary part of the complex frequency can be small enough compared with 

the Hawking temperature if the inner boundary is sufficiently close to the 

q^ I horizon, say at a proper altitude of Planck scale. Hence, the time scale of the 

instability due to the complex frequencies is much longer than the relaxation 
time scale of the thermal state with the Hawking temperature. Since ambient 
fields should settle in the thermal state in the latter time scale, the instability 
(3JT)! is not so catastrophic. Thus, the brick wall model is well defined even in a 

rotating background if the inner boundary is sufficiently close to the horizon. 
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I. INTRODUCTION 

Understanding the origin of black hole entropy is one of the most interesting problems in 
black hole physics. The black hole entropy is given by the Bekenstein-Hawking formula jl],|2j 

as 

Sbh = \A H , (1.1) 

where Ah is area of the horizon. 

It seems that full understanding of black hole entropy requires the theory of quantum 
gravity, which we do not know yet. However, we believe that general feature of black hole 
entropy can be understood by semiclassical theory, namely, quantum field theory in a fixed 
gravitational background. In fact, the brick wall model proposed by 'tHooft ||] succeeded 
to derive the proportionality of black hole entropy to the horizon area by identifying the 
black hole entropy with thermal entropy of ambient quantum fields raised to the Hawking 
temperature. It was recently clarified that in this model backreaction is small enough and 
that this model is actually a self-consistent model as a semiclassical theory |4]]. Moreover, 
it was shown that this model seeks the maximal value of entanglement entropy in the space 
of states whose backreaction is small enough |J . 

Originally, the brick wall model is proposed in the spherically symmetric, static back- 
ground, say, the Schwarzschild background. Hence, it seems interesting to see how this 
model is extended to a rotating background, say, the Kerr background |6|-§. However, it is 
known that in a rapidly rotating spacetime without horizon a field has complex-frequency 
modes [10,11 1 and that there is the so called ergoregion instability | 12fl . Thus, it might be 



expected that the brick-wall model in rotating background might be unstable and unsuitable 
for the origin of black hole entropy. 

In this paper we analyze the stability of the brick wall model in a rotating background. 
We show that the time scale of the ergoregion instability is much longer than the relaxation 
time scale of the thermal state with the Hawking temperature. In the latter time scale 
ambient fields should settle in the thermal state. Thus, the brick wall model is well defined 
even in a rotating background. 

In Sec. |I| we summarize a quantum field theory of a real scalar field in a n-dimensional 
axisymmetric stationary spacetime to show how the appearance of complex-frequency modes 
alters the structure of the quantum field theory. In Sec. |fl] we consider a scalar field in the 
4-dimensional Kerr spacetime without horizon but with an inner boundary to show the 
existence and a property of the complex frequency. Section [TV] is devoted to summarize 
this paper. 

II. SCALAR FIELD IN ROTATING BACKGROUND 

Let us consider a general n-dimensional axisymmetric stationary spacetime A4, whose 
metric is given by 

ds 2 = -N 2 dt 2 + p 2 (dip - oo B dt) 2 + q ab dx a dx b , (2.1) 



where a,b = 1, 2, • • • , {n — 2). Here, the lapse function N, the Bardeen angular velocity 
(or minus the (/9-component of the shift vector) wg, (<y2<y?)-component p 2 of the metric and 
the (n — 2)-dimensional metric q a b are assumed to depend only on the (n — 2)-dimensional 
coordinates {x a }. On this background spacetime we consider a real scalar field described 
by the action 



2 i2\ 



2J M 



d n x^)(g^d^d u <p + /xV 



(2.2) 



where the mass \x of the field can depend only on the (n — 2)-dimensional coordinates {x a }. 
We impose the boundary condition 



NfHJm^d^ = on dM, 



(2.3) 



where the boundary dM. of the spacetime M is supposed to be invariant under translations 
generated by Killing vectors d% and d^, and n M is a unit normal to dM. Note that a part of 
dM can be taken at spatial infinity. 

In this paper we quantize the system of the scalar field with respect to the time evolution 
vector D defined by 
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nM SL 



(2.4) 



where Q(x a ) is an arbitrary function of {x a } such that D is timelike in Ai. For this choice of 
the time evolution vector, it is convenient to use a new coordinate system (t, <p, x a ) defined 
by 



<£ = <£- n{x a )t 



since for this coordinate system 



D 




(2.5) 
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Following the usual quantization procedure, it turns out that the canonical momentum n 
conjugate to and, thus, the equal-time commutation relations are independent of the choice 
of tt(x a ): 
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(2.7) 



where the Lagrangian L is defined by I = J dtL. Therefore, the quantization procedure we 
follow is independent of the choice of the time evolution vector D. In this sence there is no 
ambiguity in the quantization. 

Off course, there are much freedom in selecting a ground state: we have freedom in the 
choice of a set of positive-frequency mode functions. In the following, we give one example of 
the choice of the set of positive-frequency mode functions by using a separation of variables. 
Other choices give different ground states. However, the Hilbert space of all quantum state 



is independent of the choice of the set. For example, one ground state can be expressed as 
excited states above other ground states. 

To quantize the system of the scalar field we raise the field <fi to an operator and decom- 
pose it by mode functions: 

< f>= E (®ui m aui m + Kim<ii m ) , ( 2 - 8 ) 

(ojlm)eP 



where the set P and the mode functions {<&^ m } will be defined below by Eq. (|2.16|) and 
Eq. ( |2.17| ), respectively. 

In order to define the mode functions {Q^im} in the above expansion, let us seek solutions 
{^ujim} of the field equation by the following separation of variables. 

* u im = Wx a )e- ia V m< ". (2.9) 

The function f^i m (x a ) is a solutions of the equation 
1 
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u> — Usm) rri J , 
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f uhn = 0, (2.10) 



Np^q 

with the boundary condition 

Npf ut mn"d ll f ullm = on dM, (2.11) 

where n^ denotes a unit normal to the boundary dM. (A part of dM can be taken at spatial 
infinity.) Here note that, because of the invariance of Eq. ( |2.10| ) under (u, m) <-> (—&>, —m), 
we can assume that 

JLulm = J-u*,l,—m- \^-^) 

We can choose the quantum number I so that 

J(P- 2 x^-{u;-u B m)f ulm f^ Vm = unless 1 = 1'. (2.13) 

With this choice of the quantum number I, the following property holds. 

(i&ulm, ^w'l'm^KG = unless u* = u',l = I' and m = m, (2.14) 

where the Klein-Gordon norm ($, ^)kg is given by 

($, m) KG = -i I d n - l x^-{<&Dm* - m*D§). (2.15) 

Now the set P, over which the summation is taken in Eq. (|2.8| ), is defined as 



P = P R UP C , (2.16) 

P R = {{ojlm)\u is real, (\P w j m , ^ u im)KG > 0}, 

P c = {(wIm)|Sw>0}. 



The mode functions {&ujim} in the expansion (|2.8|) are defined by 
$ ulm = — — V ulm for (ulm) G P R , 
$ ulm = — (^ m + e ia ^^^ lm ) for (ulm) G P c , (2.17) 

where the real constants C^im (> 0) and a^im are defined by 

(* w jm, ^. im )x G = C^ m e^ im . (2.18) 

For the above definition of P and {$u;«m}, the following property can be easily derived. 

(&wlm, ^w' I'm') KG = ii/^'^ram', (2-19) 

(* u l m ,& u >l>m>)KG = (2-20) 

(Cim. C'Cm')^ = —Suu,iSu'S mm i, (2-21) 

for \/(ujlm) G P and V(u/Z'm') G P. It is these properties that lead us to the above definition 
of P and {<& w im}. In Appendix [A], it is shown that the local integrability of Eq. (|2.24fj below 
requires Eq. (|2.20|) and that the normalizability of the ground state |0) requires the left 
hand side of Eq. ( [2.19| ) to be positive definite as a matrix with arguments A = (ulm) and 
A' = (oj'l'm 1 ). (Eq. (|2.21| ) is an immediate consequence of Eq. (|2.19|) .) 
Since the above orthonormality of the mode functions imply 

[a x ,a[,] = 8 xx>, [ax,a y ] = 0, (2.22) 

the Hilbert space T of all quantum states can be constructed as a symmetric Fock space 
spanned by the states |{iVx}) defined by 

\{Nx})=U { -7wr\°)^ ( 2 - 23 ) 



where the ground state |0) is defined by 

a A |0) = for VAGP. (2.24) 

Hereafter, A denotes (ulm) and A denotes (—u*,l, —m). 

The canonical Hamiltonian H [D] with respect to the time evolution vector D is given by 

H[D]= l -(D^, ( p) KG . (2.25) 

Hence, if Q is a constant Q then H [D] is a conserved quantity and can be expressed as 



1 If Q is not a constant then H[D] is not a conserved quantity in general since Qcj) does not satisfy 
the equation of motion. 



H[D] = - J2 [{&& - £lm){a x al + a[a x ) + i^u{a[a{ - a x a~ x )] . (2.26) 

Note that any states of the form (|2.23|) are not eigenstates of this Hamiltonian unless all u 
are real. (Off course, if there is no complex uj then all states of the form ( [2.23| ) are eigenstates 
of this Hamiltonian.) Moreover, in Appendix || it is shown that there is no ground state 
suitable for this Hamiltonian unless all u are real. To be precise, it is always impossible 
to eliminate terms including a\a\ or a\a\ in fl2.26|) by a Bogoliubov transformation. (See 
Appendix |TJ.) Thus, if there is a complex u then there is no stable ground state in T . 
Hence, existence of complex-frequency modes imply a kind of instability in quantum field 
theory. This conclusion is consistent with the results of Refs. [0,0 that spectrum of the 
Hamiltonian becomes continuous and that eigen states are not normalizable if there is a 
complex-frequency mode. 

Off course, there is a corresponding instability in classical theory: if there is a complex 
frequency with positive imaginary part then a solution expressed as fl2.9|) grows exponentially 
in time. On the other hand, if imaginary part of frequency is negative then the corresponding 
solution decays in time but grows exponentially in inverse time. 

Therefore, in both classical and quantum senses, appearance of complex- frequency modes 
implies instability of the system. In the next section, we show that in the brick wall model 
a scalar field has complex-frequency modes. Hence, it might be expected that the brick wall 
model might be unstable and that it might be an unsuitable model to seek entropy for an 
equilibrium state. However, it turns out that the imaginary part of the complex frequency 
can be made arbitrarily small by making the inner boundary close enough to the horizon. 
In fact, in the next section, we show that the imaginary part is small enough compared with 
the Hawking temperature if the inner boundary is sufficiently close to the horizon, say at a 
proper altitude of Planck scale. 



III. COMPLEX FREQUENCY MODES AND STABILITY OF THE BRICK WALL 

MODEL 

For simplicity, let us consider the (4 dimensional) Kerr spacetime as a background and 
suppose that the mass /x is a non-zero constant. The metric is given by 

9 / 2Mr\ , 9 S , 9 _ ,„9 „ 9 . 9 . , 9 4Mar . 9 



ds 2 = - I 1 J dt 2 + —dr 2 + Y,d6 2 + R 2 sin 2 9d<p 2 sin 2 6d<pdt, (3.1) 

where 



E = r 2 + a 2 cos 2 6, 

„2 i „2 



A = r 2 + a 1 - 2Mr, 
ZR 2 = (r 2 + a 2 ) 2 -a 2 Asm 2 6. (3.2) 



This is of the form (|2.1|) with 
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N 2 = — 
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// = iT sin z 0, 
q ab dx a dx b = ^-dr 2 + Zd6 2 . 



(3.3) 



We only consider the region r > tq in this spacetime: we impose the Dirichlet boundary 
condition on the field <fi at r = r . Hereafter we assume that r > M + -\/M 2 — a 2 : there is 
no horizon in the region r > r . 

It is well known that in this background Eq. fl2.10| ) becomes separable. In fact, we can 
find a solution of (|2.10| ) of the form 
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S = 0, 



(3.4) 



(3.5) 



and the equation for u can be written as 

d 2 u 



dx 2 



+ (y-V + )(y-V_)u = 



by introducing the non-dimensional tortoise coordinate x by 

_ x dx r 2 + a 2 
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(3.7) 



(3.8) 



Here r± = M ± V ' M 2 — a 2 , y = u;//i and V± are defined by 
2maMr 
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m 2 a 2 (r 2 + a 2 + 2Mr) 2Mr + a 2 6Ma 2 r 



(3.9) 



(r 2 + a 2 ) (r 2 + a 2 ) 3 (r 2 + a 2 ) 2 (r 2 + a 2 ) 3 

Note that in the horizon limit r —>■ r + (or x — > — oo) both of V± approach to the same value 

2\iMr + y ' 

and that V± approach to ±1, respectively, in the limit r — ► oo (or x — > oo). (See Fig. |l|.) 
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FIG. 1. The typical form of the graphs z = V±(x,y) is written on a fixed y plane. Note that 
V± depend on y as well as x through the eigen value A of equation ( |3.5| ). However, asymptotic 
behavior of V± in the limit x — > ±00 does not depend on y. 

Now let us seek complex frequency modes by examining a scattering amplitude for real- 
frequency waves. The method we shall use here is based on the following expected form of 
the scattering amplitude S near a pole y = Vr + iyi, providing that \yn\ ^> \yj\ |TT] . 
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y-VR-iyi 



(3.11] 



where 60 is a constant phase. Hence, if we can obtain this form of a scattering amplitude 
by analyzing real-frequency waves then we find an outgoing normal mode corresponding to 
V — Vr + Wi and an incoming normal mode corresponding to y = yn — iyi. After that we 
should confirm whether these normal modes converge or diverge in the limit of x — > 00. If 
these normal modes converge then they give complex-frequency mode functions. 

We first consider the case in which ma > 2/iMr + and examine the following five regime 
separately: (i) 1 < y < V_(x ); (ii) y > V+(x ); (iii) y < -1; (iv) V^(x ) < y < V + (x ); (v) 

-1 <y < i- 

(i) Ky<V-(x ) 
In this regime, let the solution of y — V- and that of y = V + be x = x\{y) and x = £2(2/), 
respectively. (See Fig. [1].) In each region separated by X\ and X2, the WKB solution is given 
by 
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for xq < x < X\, 

C 3 
ITIV4 
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/ J\T\dx) for x\ < x < x 2 , 



u 



— yjfexplij ^\T\dx)+ — -^expi-ij \j\T\dx\ for x 2 < x, (3.12) 



where 

T=(y-V+)(y-Vj. 
Hence, the standard connection formula for WKB solutions gives 
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(3.13) 
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where 
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J\T\dx 
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(3.15) 



In the limit e v — > oo, S approaches to —i unless cos£ = 0, in which case S = +i. Hence, a 
resonance will occur near a frequency corresponding to cos£ = 0. We denote the value of y 
at which 



c 



n + - I n 



by y n . Thence, we expand S near y = y n : 

y -y n + ie' 2 ' n /Aa n 
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where 
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Because of the behavior ( |3.10| ), the asymptotic behavior of a n in the limit of Xq 
be obtained 0: 
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(3.16) 



(3.17) 
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where k is 'the surface gravity of the horizon' []. This implies that a n > 0. Hence, from 
Eq. ( |3.17| ) and the last of Eq. ( p.!2| ), we can conclude that there are two regular solutions, 
whose asymptotic forms in the limit x — > oo are 



2 Although A in V± depends on y through the eigen equation (]3.5|), this asymptotic behavior of 



a n is correct since the right hand side of ( 3.10 ) is independent of A 



3 Strictly speaking, in our background there is no horizon by assumption. However, redshifted 
local acceleration is bounded from above by the surface gravity of the horizon in the extended 
spacetime which has a horizon. 
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Here the plus sign corresponds to S 1 = and the minus sign corresponds to S = 0. Thus, 
we have obtained a set of complex-frequency modes corresponding to 
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However, from the asymptotic behavior ( |3.19| ) we can conclude that 

r -r + 
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(3.22) 



where Tbh = h/2tt is the Hawking temperature of the Kerr background. 

(ii) y > V + (x ) 
In this regime, analysis depend on how many solutions y = V + has. If y = V + has no 
solution or only one degenerate solution then there is no complex frequency modes since 
whole region, x < x , is classically allowed region. \i y = V + has two solutions then we can 
repeat the above procedure for the regime (i). However, obtained WKB solutions have the 
asymptotic form ( |3.20| ) with negative a n in this case. Thus, there is no regular solutions 
which correspond to complex-frequency modes. 

(hi) y < -1 
Also in this regime, we can repeat the above procedure. If y = V- has no solution or only 
one degenerate solution then there is no complex frequency mode. For the case in which 
y = V- has two solutions, obtained WKB solution u has the asymptotic form 



u ~ exp 
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4a r 
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(3.23) 



with positive a n . Thus, there is no regular solution which corresponds to complex-frequency 
modes. 

(w)V4x )<y<V + (x ) 
In this regime, let the solution of y = V + be x = X2{y)- In each region separated by X2, the 
WKB solution is given by 
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Hence, the standard connection formula for WKB solutions gives 
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(3.25) 



(3.26) 



From this expression of S, it is evident that there is no complex-frequency mode near the 
real axis. 

(v) -1 < y < 1 
In this regime, the region with large x is classically forbidden region. Thus, there is no 
complex-frequency mode. 

Next, let us consider the case in which ma < — 2/zMr+. The above analysis can be 
applied to this case by simply replacing y with — y, V + with — V_, and V- with —V+. 
Complex frequency modes arise only in the regime V + (xq) < y < —1 and 



?B^$fw ~ ±ne~ 2v 



in r°' r+ 



(r ^r + ). (3.27) 



Finally, let us consider the case in which —2[iMr + < ma < 2/iMr + . From the above 
analysis for other cases, it is evident that there arise no complex-frequency mode functions 
since neither the regime 1 < y < V-(xo) nor the regime V + (xq) < y < —1 exist in this case. 

In summary, in this section we have shown that in the Kerr background without horizon 
but with an inner boundary a scalar field has complex-frequency modes and that the imagi- 
nary part of the complex frequency is small enough compared with the Hawking temperature 
if the inner boundary is sufficiently close to the horizon, say at a proper altitude of Planck 
scale. 

IV. SUMMARY AND DISCUSSION 

We had analyzed the stability of the brick wall model in a rotating background. We 
had shown that in the Kerr background without horizon but with an inner boundary a 
scalar field has complex- frequency modes and that, however, the imaginary part of the 
complex frequency can be small enough compared with the Hawking temperature if the inner 
boundary is sufficiently close to the horizon, say at a proper altitude of Planck scale. Hence, 
the time scale of the ergoregion instability is much longer than the relaxation time scale of 
the thermal state with the Hawking temperature. In the latter time scale ambient fields 
should settle in the thermal state. In this sense ergoregion instability is not so catastrophic. 
Thus, the brick wall model is well defined if the inner boundary is sufficiently close to the 
horizon. 

Now, let us discuss physical interpretation of the existence of complex-frequency modes. 
First, for a rotating black hole background, there is no complex frequency mode |TJ]]. How- 
ever, it is well known that superradiant modes of fields are amplified by scattering. On the 
other hand, for the brick wall (or a rapidly rotating star) background, the amplification of 
superradiant modes can be suppressed by a boundary condition say, the Dirichlet boundary 
condition at the inner boundary [IT)]. In stead of the amplification of superradiant modes, 



as shown in this paper, for this background there appear complex frequency modes. These 
complex frequency modes are outgoing and incoming normal modes of the field and may be 
understood intuitively as "quasi-bound states" in the ergoregion. (See Eq. fl3.16|) .) This in- 
terpretation is consistent with the fact that for a black hole background there is no complex 
frequency mode since there is no " quasi-bound state" in the ergoregion: any excitations with 
negative energy w.r.t. observers at infinity will fall into the hole. Based on this observation, 

11 



thus, it is expected that in the brick wall background the imaginary part of the complex 
frequency should become arbitrarily small in the limit that the inner boundary becomes 
close enough to the horizon since in this limit there appears a large room for the excitations 
with negative energy w.r.t. observers at infinity to escape to. This consideration is, off 
course, consistent with our result in this paper: we have shown that the imaginary part of 
the complex frequency can be small enough compared with the time scale determined by 
the Hawking temperature if the inner boundary is sufficiently close to the horizon, say at a 
proper altitude of Planck scale. 

Next, let us discuss a relation to the so called Schiff-Snyder- Weinberg effect [IB]. In 



Ref. |T7| the relation between the Klein paradox []T8j] and superradiance in a rotating black 
hole background was discussed in detail by using a rectilinear model of the Kerr spacetime. 
Since the situation in the Klein paradox can be understood as a limit (the so called Klein 
limit) of the Schiff-Snyder- Weinberg effect [17], the situation in a rotating black hole back- 



ground should be understood as a limit of our situation, i.e. the brick wall model. This is 
actually the case: in the Xq — > — oo limit the complex frequencies in the brick wall back- 
ground disappear as in the Kerr black hole or in the rectilinear spacetime considered in 
Ref. 0. 

Finally, we would like to mention a possibility to stabilize the quantum field theory in a 
rotating background by introducing a nonlinear interaction. In Refs. [B],2C] it was suggested 



that a nonlinear interaction will prevent the vacuum from being unstable even if there are 
complex-frequency modes. It will be interesting to investigate such a possibility in the case 
of quantum field theory in a rotating background. Physics in a rapidly rotating background 
spacetime will be as interesting as physics of strong fields |21| . 
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APPENDIX A: KLEIN-GORDON NORM AND INTEGRABILITY 

In Sec. |IT] we have expanded the field operator as Q2.8p by mode functions {$ui m } 
satisfying (|2.19| - |2~2~T|) . In other words, we have required that coefficients of annihilation 



operators should have positive Klein-Gordon norm instead of requiring positivity of the 
frequency u> with respect to the Killing time t. 

In this appendix we show that the positivity of the Klein-Gordon norm is required by 
integrability of equations for the ground state, say, Eq. ( |2.24| ). 



Let us consider a scalar field described by the action (2.2) in a general n-dimensional 
globally-hyperbolic spacetime: 

ds 2 = -N 2 dt 2 + 1]k (dx J + (3 3 dt){dx k + (3 k dt), (Al) 
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where the lapse function N, the shift vector /3 J , the (n — l)-dimensional metric 'jjk and the 
mass \x can depend on both t and the (n — l)-dimensional coordinates {x- 7 }. 

To make our arguments definite, let us discretize the system of the scalar field. Since 
the Lagrangian L defined by I = J dtL can be written as 



l = ]- I <r- v r 



^(d t <p - pd^y - N^i^d^d^ + /i 2 2 ) 



(A2) 



we can discretize it to obtain f] 



L = \G AB {^ A - f£<P c )(<P B - fE<P D ) - \V AB <P A <P B , (A3) 

provided that we suppose the following correspondence. 

' (F^xN^fi^drfd^ + /x 2 2 ) <s> V AB (f> A 4> B . (A4) 



These relations will be used when we take a continuous limit. The corresponding equation 
of motion and the commutation relations are 

(G AB <j> B - G AB fE<t> C ) + G B cfy B + (Vab ~ f A G C DfE)^ B = 0, (A5) 

and 

[0 A , n B ] = i6 A , [c/> A , <f> B ] = [n A , n B ] = 0, (A6) 

where the momentum n A conjugate to (fi A is defined by 

BT 

7r A = ^ = G AB (<p B -fS<P c ). (A7) 

Now let us expand the field operator by "mode-functions" {§ A , & A *} (n — 1, 2, • • •), i.e. 
solutions of the equation of motion: 

^ = EMn+^f), (A8) 

n 

where we assume that {§ A , § A } forms a complete set of linearly independent solutions of 
the equation of motion. 

We would like to define the corresponding ground state by 



4 Although explicit forms of the matrices G, f and V depend on the way of discretization, all 



we need to take the continuous limit in the following arguments is the correspondence (A4) only. 
Thus, the result in this appendix is independent of an explicit way of discretization. 
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o n |0) = 0, for Vn. 



(A9) 



This is the discretized version of Eq. (12.24]). However, this equation is not integrable in 



general. Hence, in this appendix we would like to seek the necessary and sufficient condition 
for the integrability of this equation. 

First, it is easily shown that a certain linear combination of (p A and tta is written as a 
linear combination of a„ as follows. 



■Ka — i&>AB< 



lB 



J2an$n(tt + tt*)AB, 



(MO) 



where the matrix Qab is defined by 

Thus, in order for Eq. ( |A"S| ) to be integrable, it is necessary that 

l(n A - in AC <l> c ), (ti b ~ «£W D )]|0> = {^ab ~ n BA )\0) = 0, 
which is equivalent to 



(AH) 



(A12) 



$f (Pab - <W« = -iGAB[$f(K* - fE*%) - (K* - fcK*)Kl = 0. (A13) 



Next, by using the relation ( |A10| ), the solution of Eq. ( |A9| ) is obtained at least locally as 
follows, provided that the local integrability condition (|A13| ) is satisfied. 



({^}|0)=ATexp 



— ^Qab< 



kAiB 



(A14) 



In order for this wave function to be well-defined, i.e. normalizable, it is necessary and 
sufficient that the hermite matrix (Q + Q^ab be positive definite. This condition is restated 
that the matrix X nm defined as follows should be positive definite. 



x nm = ^(n + rt) AB ^ = -iGAB&H** - /£*£•) - iK - f£*Z)v, 



(A15) 

In summary, in order for Eq. ( |A9[ ) to be integrable, it is necessary and sufficient that the 
condition ( |A13|) is satisfied and that the matrix X nm defined by (|A15 ) is positive definite. 
These two conditions can be restated as follows: 



where the norm ($, \I/) is defined by 



0, for Vn, m, 
is positive definite, 



($, *) = -iG A b[® A {^ B * ~ fE^ C * 



($ A -/£$ G )# 



r-B*l 



(A16) 

(A17) 



(A18) 



) that this norm is constant in time 



It is easy to show by using the equation of motion 
if both $ and \l/ satisfy the equation of motion. 

Provided that the condition ( A17|) is satisfied, it is possible to take linear transformation 
of {& n } so that 
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($„,$ m ) = 5 nm , (A19) 

preserving the condition (|A16|). In this normalization, it can be easily shown that 



|a n , a m J — o nm , [a n , a m \ — [a n , a m \ — (J. (A2UJ 

Thus, we can construct the Hilbert space of all quantum states spanned by 

IW>=(n%-J|0>- (A21) 



Now let us take a continuous limit. From the correspondence flAJ), it is evident that the 
norm (|A18|) is a discretized version of the Klein-Gordon norm 



($, V) KG = -iJ d^x^/ju^d^* - #*<^$), (A22) 

which reduces to ( |2.15|) for a spacetime metric of the form (|2.1|) . 

It is also evident that the integrability conditions ( |A16| ) and ( |A19| ) in the continuous 
limit for the expansion (|J) is that (^20|) and ( gig ) for \/(ulm) G P and \/(uj'l'm') G P. 

APPENDIX B: HAMILTONIAN FOR A COMPLEX FREQUENCY MODE 

In this appendix, we show that it is always impossible to eliminate terms including a\a\ 
or a\a\ in ( |2.26| ) by a Bogoliubov transformation. 

In general, a Bogoliubov transformation can be written in terms of two matrices a and 
j3 satisfying 

QJQjt _ = ! 

af - l3a T = (Bl) 

as 

$ n -> ^(«nm$m + Pnm®* m ), ( B2 ) 

m 

where {$ n } {n — 1, 2, • • •) is a set of positive frequency mode functions. Let us consider a 
Hamiltonian of the form 

H = ~E nm (a n al + 4a m ) + ^A nm a n a m + ^A nm *al al, (B3) 

where E is a hermite matrix and A is a symmetric matrix. Under the Bogoliubov transfor- 
mation, the coefficient-matrices E and A are transformed as 



V ; 



E -* (a£a f + /3£*/3 f ) + (aA/3 f - /?Aa 

A -> (aAa T - (3A*f3 T ) + (aEp T + f3E*a T ) . (B4) 
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Returning to the problem, the contribution of a pair of complex-frequency modes A and 
A to the Hamiltonian (|2.26|) can be written as 



h = -(9?co> — Vim) (a\a' x + a x a\) — (a x a' x + a' x a x ) + i^suj(a x a' x — a\a x ) 
This is of the form (IB 31) with 



(B5) 



E = (fcu - fim) 



1 
-1 



A = — i^u 



1 

1 



(B6) 



Hence, what we shall show now is that there is no choice of 2 x 2 matrices a and f3 satisfying 
Eq. (|BTD and 



(aKa T - pA*p T ) + (aE(3 T + pE*a T ) = 0. 



(B7) 



This statement is easy to show. First, since the first of Eq. ( |B~T| ) implies that a has the 
inverse, the second of Eq. (|B1| ) and Eq. ( |B7|) are written as 7 = ■y T and 



iQuj 



1 

1 



AVoh T 



($tu - Vim) 



1 ° 7 r +7 r ° 

-1 ) 1 '0-1 



(B8) 



where 7 = a~ l (3. These equation are easy to solve with respect to 7. The result gives 
I det 'y I = 1. However, from the first of Eq. ( |B~T| ) it is derived that | det -y | < 1, which 
contradicts with the above result. Therefore, there is no choice of 2 x 2 matrices a and (3 
satisfying Eq. QBTp and (§7|). 
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